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Abstract: 

Crucial information on the coupling of the Higgs boson to bottom quarks is expected from 
Higgs production in association with a forward tagging jet at a future high-energy lepton-proton 
collider. In order to control the theoretical uncertainties of the signal process, the impact of 
radiative corrections has to be quantified. We present the full next-to-leading order QCD correc- 
tions to e~p — > e~jH and e~p — > v e jH in the form of a flexible Monte-Carlo program allowing 
for the calculation of cross sections and kinematic distributions within experimentally feasible 
selection cuts. QCD corrections are found to be very small for cross sections, while the shape 
distortion of distributions can be as large as 20%. Residual scale uncertainties at next-to- leading 
order are at the permille level. 



1 Introduction 



With the start-up of the CERN Large Hadron Collider (LHC), promising prospects have arisen 
for pinning down the mechanism responsible for electroweak symmetry breaking. In order to 
verify the scenario realized by the Standard Model (SM), a neutral, CP-even Higgs boson has 
to be discovered and its properties have to be investigated carefully. This implies measurement 
of its couplings to gauge bosons and fermions. At the LHC, the determination of the Higgs 
couplings to the top quark, tau lepton, and the weak gauge bosons should be possible, once a 
sufficient amount of data has been collected [IHS]. If the Higgs boson turns out to be relatively 
light, its dominant decay mode is to bottom quarks. Extracting the bottom Yukawa coupling 
is of particular importance, as this quantity will help to distinguish the SM from various of its 
extensions such as supersymmetric scenarios. 

However, probing the Hbb coupling at the LHC is a challenging enterprise due to large QCD 
backgrounds in the dominant production modes gg — > H, associated ZH and WH production, 
and tiH production with subsequent decay H — ?■ bb Higgs production via weak-boson 
fusion (WBF) with decay into bottom quarks, on the other hand, provides challenges for the 
ATLAS and CMS triggers [I]. Therefore, additional search channels, such as WBF production of 
7-ff [5] and WH [6], with H — > bb have been considered. Alternatively, one could make use 
of the substructure of so-called "fat" jets resulting from the decay products of heavy particles 
in different production modes [THS]. Despite these recent developments, the observation of the 
H — > bb decay mode remains challenging at the LHC. 

New opportunities could be provided by a future high-energy lepton-hadron collider such as 
the CERN Large Hadron electron Collider (LHeC), which is currently under scrutiny [10] . For 
now, various scenarios are being discussed, which all make use of the LHC proton beam with an 
energy of E p = 7 TeV and an electron beam in the range of E e = 50 — 200 GeV, corresponding 
to a center-of-mass (c.m.s.) energy of \fS = 2^jE p E e w 1.18 — 2.37 TeV. Several studies have 
been performed to estimate the capability of such a machine to access the coupling of the Higgs 
boson to bottom quarks in a cleaner environment than at a hadron-hadron collider [11]|12]. In 
particular, Han and Mellado [13J provided a detailed, cut-based leading-order (LO) analysis of 
the Higgs signal at the LHeC in the presence of various backgrounds. 

In addition to identifying the selection criteria most suitable for an extraction of the Higgs 
signal from backgrounds at an ep collider, precise knowledge of the signal processes itself is 
crucial. We therefore provide accurate predictions for cross sections and kinematic distributions 
of Higgs production in association with a tagging jet at the LHeC, taking next-to-leading order 
(NLO) QCD corrections fully into account. We start with a brief outline of the calculational 
methods used in Sec. [2j In Sees. 13.11 and 13.21 we present our numerical results for the charged 
current and neutral current production modes, respectively. Our focus is the estimate of the 
perturbative uncertainties associated with predictions for the signal process. Conclusions are 
given in Sec. HI 

2 Setup of the calculation 

At LO, Higgs production in ep collisions mainly proceeds via lepton-quark scattering, mediated 
by the exchange of massive weak gauge bosons which in turn radiate off a Higgs boson. In 
the final state, the scattered quark gives rise to a jet j. Depending on the charge of the weak 
exchange boson, the respective production modes, ep — > v e jH and ep — > ejH, are referred to 
as charged current (CC) and neutral current (NC), respectively. For simplicity, we will focus on 
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Figure 1: Leading order Feynman diagram for Higgs production at a lepton-hadron collider in 
the neutral current mode eq — > eqH. 



the NC mode for the discussion of the technical details of the calculation. The treatment of the 
CC processes proceeds on the very same lines. 

Figure [1] depicts the single Feynman diagram for the dominant NC production mode at parton 
level, 

e(h)q(p b ) -> e(k 2 )q( Pl )H(P H ). (2.1) 

The anti-quark initiated process, eq — > eqH, is easily obtained thereof by crossing relations. The 
scattering amplitudes, -M^q and A^q, for each partonic subprocess are computed numerically 
with the helicity-amplitude techniques of [14015] . For the hadronic jet cross section, each partonic 
matrix element squared, |.Mlq| 2 , nas to be convoluted with the parton distribution function of 
parton b in the proton, fbixb,^), multiplied with the flux factor and a suitable jet-defining 
function, and integrated over the available phase space, 

a LO (ep^ ejH)= V f dx b f q (x b , Mf ) 1 d$ 3 \M b LO \ 2 F^ fa) . (2.2) 

Here, Xb is the fraction of the proton momentum taken by parton b at a factorization scale fi-p. 
The partonic c.m.s. energy squared is given by s = {k\ + Pb) 2 , d$> m denotes the m-particle phase 

(n) 

space of the system, and the function F) defines the jet algorithm for an n-parton final state. 

At NLO, virtual and real-emission corrections to the tree- level amplitudes have to be con- 
sidered. Since at the perturbative order we consider the color-neutral lepton-beam particles are 
not affected by QCD corrections, the structure of the radiative corrections is simple and resem- 
bles that for deeply inelastic lepton-nucleon scattering. Very similar contributions arise in H jj 
production via WBF in hadron-hadron collisions and have been computed in Refs. |16H18j . For 
our work, we adapt the NLO-QCD corrections of Ref. [16J to Higgs production at lepton-hadron 
colliders. 

The real-emission amplitudes, M h R: are obtained from the tree-level diagrams by attaching 
one gluon to the quark line in all possible ways. The resulting diagrams are depicted for the 
representative quark-initiated subprocess, 

e(h) q( Pb ) e(k 2 ) qfa) g(p 2 ) H(P H ) , (2.3) 

in Fig. [21 In addition to the eq and eq channels with an extra gluon in the final state, now also 
the gluon-initiated subprocess, 

e(h)g( Pb ) e(k 2 ) qfa) q{p 2 ) H(P H ) , (2.4) 
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Figure 2: Real emission contributions to Higgs production at a lepton-hadron collider in the 
neutral current mode eq — > eqgH. 
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Figure 3: Real emission contributions to Higgs production at a lepton-hadron collider in the 
neutral current mode eg — > eqqH. 



has to be considered, cf. Fig. For each channel, in regions of phase space where soft and/or 
collinear configurations can occur, singularities are encountered in the phase-space integrals of the 
real-emission amplitudes squared. These divergences are regularized dimensionally in d = 4 — 2e 
space-time dimensions. The cancelation of the singularities is performed in terms of the dipole 
subtraction formalism of Ref. [TU]. With the radiative corrections being of the same form as 
in WBF at hadron colliders, the relevant subtraction terms, |-A4g ing | 2 , can be adapted from 
Ref. |16| by replacing the entering Born-type parton-parton amplitudes with the appropriate 
lepton-parton contributions. The finite 2-parton cross section for the NC eq subprocess is then 
of the form 



a. 



NLO,<? 



dx b f q (x b , fiy) — d®4 
is 



\M q R \ 2 Ff\ Pl , P2 )-\M 



q 

sing 



(2.5) 



where, in the notation of [19] , p\ corresponds to a set of parton momenta approaching Born 
kinematics in the collinear and soft limits, pi — > p\. In these regions, the jet-defining function 



has to obey F 



(2) 



,(1) 



Expressions similar to Eq. (|2.5p are obtained for the anti-quark- and 



the gluon-initiated subprocesses. 

In order to arrive at well-defined results for physical observables, the counter-terms that have 
been introduced in Eq. (|2.5p are integrated analytically in d dimensions over the phase-space of 
the potentially soft or collinear parton, yielding 



(1(e)) = \M 



q 1 2 <Xs(PR) n 

Lo1 ~^T Cf 



47T/4 

Q 2 



r(l + e) 



+ - + 9 - -TV 2 



(2.6) 



where a s (//R,) is the strong coupling as a function of the renormalization scale /xr, Cf = 4/3, and 
Q 2 = —(pb — Pi — P2) 2 denotes the virtuality of the weak boson attached to the quark line. The 
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Figure 4: Virtual contribution to Higgs production at a lepton-hadron collider in the neutral 
current mode eq — > eqH. 



virtual corrections comprise vertex corrections to the quark line as depicted in Fig. HJ They are 
entirely proportional to the Born amplitude such that the relevant interference with the tree-level 
amplitude is given by 



2Re[Ml / MlQ = \M[ \ 2 ^^- 



C F 



4"7T/i 



Q 2 



r(i + £ ) 



e 2 e 



+ c 



virt 



(2.7) 



where c v i rt stands for a constant given by c v i r t = n 2 /3—7 in dimensional reduction (c v i r t = 7r 2 /3— 8 
in dimensional regularization) . After adding the integrated counter-terms, the resulting 1-parton 
contribution to the NLO cross section is finite and for the NC eq subprocess takes the form 



a 



NLO,g 
1 



jf 1 dx b f g (x b , MF ) 1 d<S> 3 \Ml Q \ 2 Fj V) 1 + C F U - -^ 2 



+ c- 



virt 



.(2. 



A remaining divergent piece of the integrated counter terms is factorized into the distribu- 
tion function of the incoming parton. For the quark-initiated contribution, the surviving finite 
collinear term is given by 



a 



NLO,<? 
lxoll 



f 1 dx b f c q (x b ,fi F ,fi R ) id$3 \M q LO \ 2 F^\ Pl ) , 
Jo Zs 



(2.9) 



where f^(x b , /j, F , I^r) is defined by 



fC ( \ «s(mr) f 1 dz t fx \ 



+ fq - Zfq {%,I*f) B(z) + f q (^-,H F ^JC(z)} 

o s (^r) 



+ 



2tt 



-fq{x,fJ>F)D(x) 



with the integration kernels 

A(z] 



T F [z 2 + (1 - z) 2 ] In Q2{ ~\ Z) + 27) : ( 1 : ) 



,B(z) = C F 
C{z) = C F 
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Q 2 



2 4{l 



Q 2 

/4 



27 

+ 21n(l-x) ln^2-+ln 2 (l-x) + vr 2 - — - c virl 



(2.10) 

(2.11) 
(2.12) 
(2.13) 
(2.14) 
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where Tp = 1/2. Analogous expressions are encountered for the anti-quark initiated subpro- 
cesses. 

The NC cross-section contributions discussed above and the respective CC expressions have 
been implemented in a flexible Monte-Carlo program that allows for the calculation of experi- 
mentally accessible observables within realistic acceptance cuts at NLO-QCD accuracy. In order 
to ensure the reliability of our calculation, a number of tests has been performed: 

• The tree-level and the real-emission amplitudes have been compared to the corresponding 
expressions generated automatically by MadGraph [20] for a representative set of phase- 
space points. We found full agreement. Since the virtual amplitudes are just multiples of 
the Born cross section, they do not need to be tested separately. 

• The integrated LO cross sections have been compared to the corresponding results of the 
MadEvent package [Hl[22] for the inclusive selection cuts to be described in Sec. [3j The 
cross sections agree within the numerical accuracy of the two programs. 

• The real-emission matrix elements have been found to vanish when the polarization vector 
of the external gluon is replaced with its momentum. This procedure tests the QCD gauge 
invariance of the real-emission amplitudes. 

• To probe our implementation of the dipole subtraction, we checked that the real-emission 
contributions approach the subtraction terms in singular regions of phase space. 

3 Numerical results 

In order to assess the impact of NLO-QCD corrections on cross sections and kinematic distribu- 
tions, we consider Higgs production at the LHeC within different settings. We present numerical 
results applying only minimal selection cuts first and then turn to the settings which have been 
suggested in Ref. [13] for extracting the Higgs signal from various backgrounds at the LHeC. 
Unless stated otherwise, we assume electron-proton collisions with 



and a SM-like Higgs boson with a mass of Mh = 120 GeV. Quark and lepton masses are set to 
zero throughout, and contributions from incoming top and bottom quarks are neglected. 

To simulate a generic Higgs decay without specifying a particular channel, we generate an 
isotropic Higgs boson decay into two massless particles (which represent 77 or bb final states), 
requiring each decay particle, labeled d, to be separated from a jet in the rapidity-azimuthal 
angle plane by 



The respective branching ratio BR{H — > dd) is not included in the numerical results we present 
below. 

Depending on the order in perturbation theory at which an observable is being computed, 
we use the LO or the NLO set of the MSTW08 parton distribution functions (PDFs) [23]. For 
the strong coupling, which enters only at NLO, we take the appropriate two- loop expression. 
We have chosen the weak gauge boson masses, Mw = 80.423 GeV, Mz = 91.188 GeV, and the 
Fermi constant, Gp = 1.166 x 10 -5 / GeV 2 , as electroweak input parameters. Thereof, oqed 
and sin 2 9\y are computed via LO electroweak relations. Final-state partons are recombined into 



E e = 140 GeV and E p = 7 TeV , 




AR jd > 0.4 . 
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Figure 5: Dependence of the ep — > v e jH cross section, Cqqj at the LHeC on the factorization 
and renormalization scales. The NLO curves show (Tqq clS cl function of the scale parameter £ 
for /in = p-p = £p (solid red). The LO cross section only depends on pp = t;p (dashed black). 



jets according to the algorithm [24-27J with a resolution parameter D = 0.4. Jets are required 
to have a transverse momentum of 

p£* > 15 GeV , (3.3) 

and events which do not exhibit at least one jet are disregarded. In configurations with more 
than one identified jet, the jet of highest transverse momentum is referred to as "tagging jet". 
Below, we will refer to observables within the cuts of Eqs. (|3,2p and (|3.3p as "inclusive" . 
For clarity, we will discuss the CC and NC production modes separately in the following. 

3.1 Charged current processes 

CC processes represent the major source of Higgs bosons in ep collisions due to the value of the 
charged weak boson coupling to the electron, which is much larger than the respective neutral 
current coupling. At LO, the signal topology consists of large missing energy from the neutrino, 
a dd pair from the Higgs decay, and a forward energetic tagging jet. At NLO, this jet may be 
composed of two partons, or we may encounter a second, well-separated jet. 

In Fig. [5j the integrated CC cross section within the cuts specified in Eqs. (|3.2p and ()3.3|) . 
ctqq, is displayed as a function of the scale parameter £, which is related to the factorization and 
renormalization scales, pp and /Ur, via 

Mf = £ Mo and ^r = £ p with p = Q . (3.4) 

The LO cross section depends only on pp, and the variation of (Tqq with £ resembles the scale 
dependence of the quark distribution function f g (x, /if) at relatively large values of x. At NLO, 
the //F-dependence of the parton distributions is partly canceled by respective contributions in the 
partonic cross section. The renormalization scale enters via the strong coupling a s (p,^). Due to 
the small size of the NLO-QCD corrections, the overall dependence on jUr is small. Qualitatively, 
the improvement of the scale dependence when going from LO to NLO is very similar to what 
has been observed for WBF processes in hadronic collisions [28H3T] . It is interesting to note, 
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however, that the relative size of the scale uncertainties is much smaller in the ep case, due to 
the presence of only one parton distribution function as opposed to two in pp collisions. 

Additional uncertainties could arise from the incomplete knowledge of the parton distribution 
functions of the proton. In order to estimate this uncertainty, we have calculated (Tqq for two 
scale choices, Hf = HR = Q and = PR = Mw, using two different PDF parameterizations: 
our default set MSTW08 and the CTEQ6 set [32]. In particular, the CTEQ6M set is used at 
NLO, and the CTEQ6L1 set at LO. Our results are listed in Table [TJ Very similar results are 



PDF set 


^cc' LU (/«o = M w ) 


inc,lNLO/ tit \ 
a CC CA*0 = M W) 


inc,LO/ /~,\ 

^cc Wo = Q) 


inc,JNLO/ 

°"cc w = Q) 


CTEQ6 


241.83 fb 


240.82 fb 


241.89 fb 


239.78 fb 


MSTW08 


240.11 fb 


241.35 fb 


240.24 fb 


240.30 fb 



Table 1: Inclusive cross section, (Jqq, for ep —> u e jH at the LHeC for the scale choices fio = Mw 
and hq = Q with £ = 1, and for two different sets of parton distribution functions. The statistical 
errors of the quoted results are at the sub-permille level and therefore not given explicitly. 

obtained when Mw is used as a default scale rather than Q, since Q tends to be close to the weak 
boson mass scale in this class of reactions. The numerical uncertainty of (Jqq due to the PDF 
parameterization is small both at LO and NLO, with relative differences between the two choices 
being at the permille level. In the following we will stick to the MSTW08 parton distributions 
and set fj,F = Mr = Q, unless stated otherwise. 

In order to estimate the impact of the NLO-QCD corrections on various observables, we 
introduce the dynamical iT-factor of a distribution da/dO, 

_ da^°(HK,HF)/dO 
K[X) ~ do^M/dO ■ 

In Fig. [6l o"q§ is shown as a function of the Higgs boson mass together with the associated K- 
factor for two different choices of the scale parameter £ with no = Q- The NLO-QCD corrections 
are small over the entire mass range considered with relative scale uncertainties increasing with 
Mff. Fig. [7] illustrates that the scale uncertainties go down when the energy of the electron beam 
is increased. Considering this feature together with results on the signal significance reported in 
Ref. [13], we conclude that for a given Higgs mass large values of E e are favorable for a clean 
extraction of the Higgs signal at the LHeC. In the following, we will stick to our default choice, 
E e = 140 GeV and M H = 120 GeV. 

Having estimated theoretical uncertainties for ctqq we now turn to a study of CC Higgs 
production at the LHeC with selection cuts for an optimal signal significance, following Ref. [13] . 
For each decay product d of the Higgs boson we require 

p d T > 30 GeV , \r] d \ < 2.5 , AR jd > 0.4 , (3.6) 

with pj, denoting its transverse momentum and rj d its pseudo-rapidity being given in the labo- 
ratory frame of the ep system. The neutrino gives rise to missing energy in the detector, which 
we request to obey 

Ef ss > 25 GeV . (3.7) 

Any identified jet has to fulfill 

pt* > 15 GeV and |y jet | < 5, (3.8) 
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Figure 6: Panel (a) shows the dependence of the ep — > v e jH cross section, (Tqq, at the LHeC 
on the mass of the Higgs boson at LO (black dashed) and NLO (red solid) for E e = 140 GeV. 
In panel (b) the corresponding .ftT-factor as defined in Eq. ()3.5[) is shown for different values of 
the factorization and renormalization scales, /2r = //p = £ Q- 




Figure 7: Panel (a) shows the dependence of the ep —> v e jH cross section, c^qi at the LHeC on 
the energy of the electron beam at LO (black dashed) and NLO (red solid) for Mh = 120 GeV. 
In panel (b) the corresponding i^T-factor as defined in Eq. ()3.5[) is shown for different values of 
the factorization and renormalization scales, /^r, = fip = £Q. 
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Figure 8: Transverse momentum distribution of the tagging jet in ep — > v e jH at the LHeC at 
LO (dashed black) and NLO (solid red) [panel (a)] and associated i^-factor [panel (b)] as defined 
in Eq. ()3.5[) for different values of the factorization and renormalization scales, /iR = fj,p = £Q- 



where y jet denotes the rapidity of a jet. For the tagging jet we furthermore demand 

p T ag > 30 GeV and - 5 < y tag < -1 . (3.9) 

Note that (contrary to Ref. [13]) we count positive rapidity in the direction of the electron beam. 
The invariant mass of the Higgs boson candidate and the tagging jet, Mn.tag, is required to be 
large, 

M H ,ta g > 250 GeV . (3.10) 

The CC cross section within the cuts given in Eqs. (13.6p - (13.10j) . crg^, changes by about 6.5% at 
LO and 8%o at NLO when //r and //p are varied simultaneously in the range Q/2 to 2Q. The scale 
dependence of the CC cross section is thus slightly larger after the application of selection cuts 
than in the inclusive case, but still very small. In particular, the NLO prediction is extremely 
stable with respect to scale variations irrespective of the selection cuts applied, which indicates 
that Higgs production at the LHeC is under excellent control perturbatively. 

Similarly to WBF processes at hadron colliders, Higgs production processes at the LHeC 
feature tagging jets with pronounced kinematic properties that help to extract the signal from a 
variety of backgrounds. In ep — > v e jH, the tagging jet tends to be located in the far-backward 
region of the detector, along the direction of the incoming proton, and can be well separated 
from the decay products of the Higgs boson. This is particularly important, if Higgs production 
at the LHeC is to be utilized for an extraction of the coupling of the Higgs boson to bottom 
quarks. In this context, it is important to estimate the impact of NLO-QCD corrections on 
kinematic distributions of the tagging jet, such as its transverse momentum distribution displayed 
in Fig. [8] for fJ-R = fJ-F = Q together with the associated if-factor for two different values of the 
factorization and renormalization scales. The corrections to da/dp^ g are slightly positive at low 
Py g , but increasingly negative towards larger transverse momenta. The difference between the 
-fT-factors for the two considered values of £ indicates the scale uncertainty of the distribution, 
which also increases with p^ g . 
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Figure 9: Dependence of the ep — > ejH cross section, cx^c; at the LHeC on the factorization 
and renormalization scales. The NLO curves show function of the scale parameter £ 

for hyi = = £mo (solid red). The LO cross section only depends on /xp = £(j,q (dashed black). 



3.2 Neutral current processes 

Cross sections for NC Higgs production via ep —> ejH are considerably smaller than the related 
CC quantities. Due to the charged lepton in the final state, however, the event reconstruction 
in the NC mode is superior to the CC processes with a neutrino that can only be identified 
indirectly as missing energy in the detector. Moreover, NC processes yield information on the 
coupling of the Higgs boson to neutral gauge bosons rather than charged ones which are probed 
in CC Higgs production at the LHeC. 

The theoretical uncertainties of the ep — > ejH cross section at the LHeC are very similar 
to those encountered in the CC production mode. Figure [9] shows the scale dependence of the 
NC cross section, c^C' at the LHeC within the cuts of Eqs. (13. 2 j) and (|3.3j) as function of the 
scale parameter £, defined by Eq. (|3.4p . As before, we use no = Q a default scale. In Fig. [TOl and 
Fig. [TTT c^c is shown as a function of the Higgs boson mass and of the energy of the electron 
beam, respectively. The corresponding K-factors illustrate the associated scale uncertainties 
when and /ip are varied simultaneously in the range Q/2 to 2Q. Similarly to the CC case, 
the NLO-QCD corrections are small, but increase with Mh, while they decrease when E e goes 
up. 

In Ref. [13], selection cuts have been presented for separating the NC Higgs signal from 
the leading background processes, in particular e~p — > e~bbj + X. We apply the same set of 
cuts in addition to the inclusive cuts of Eqs. (13.21) and (|3.3h . assuming Mh = 120 GeV and 
E e = 140 GeV. Any jet has to be located within the rapidity range 

|y jet |<5. (3.11) 

Furthermore, the tagging jet has to exhibit a transverse momentum larger than 

Py g >30GeV. (3.12) 

The scattered electron is supposed to obey 

p e T > 30 GeV and |r/ c | < 5 , (3.13) 
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Figure 10: Panel (a) shows the dependence of the ep — > ejH cross section, cr^g, at the LHeC 
on the mass of the Higgs boson at LO (black dashed) and NLO (red solid) for E e = 140 GeV. 
In panel (b) the corresponding ET-factor as defined in Eq. (13.5 [) is shown for different values of 
the factorization and renormalization scales, fJ.R = fJ-F = C Q ■ 




Figure 11: Panel (a) shows the dependence of the ep — > ejH cross section, cr^g, at the LHeC on 
the energy of the electron beam at LO (black dashed) and NLO (red solid) for Mh = 120 GeV. 
In panel (b) the corresponding ET-factor as defined in Eq. (I3.5P is shown for different values of 
the factorization and renormalization scales, /^r, = = £,Q- 
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Figure 12: Transverse momentum distribution of the scattered electron in ep — > ejH at the 
LHeC at LO (dashed black) and NLO (solid red) [panel (a)] and associated if-factor [panel (b)] 
as defined in Eq. (I3.5P for different values of the factorization and renormalization scales, /jr, = 



while for the decay products of the Higgs boson we require 

p% > 30 GeV and \r] d \ < 2.5 . (3.14) 

With this set of selection cuts, the NLO-QCD corrections to the ep — > ejH cross section, a^f, 
are very small, amounting to about 2% for //r = /xp = Q. At NLO, cr™^ changes by only about 
2%o when the scales are varied in the range Q/2 to 2Q. Nonetheless, the shapes of kinematic 
distributions can receive larger corrections, as we will see below. 

We first explore the kinematic properties of the scattered electron. Figure H2] (a) shows 
the transverse-momentum distribution of the electron in the final state. The shape of this 
distribution is barely affected by NLO-QCD corrections. The corresponding X-factor is shown 
for /xf = fJ>R = Q/2 and 2Q in Fig. [12] (b), with the difference between the curves indicating 
the scale uncertainty of the prediction. The relative scale dependence is small, but increases 
with pj,. Similarly to da/dp^, the rapidity distribution of the electron retains its form at NLO. 
Radiative corrections do have an impact on jet observables, however. In Fig. [13] the rapidity 
distribution of the jet with the highest value of |y jet | is shown together with the corresponding 
-RT-factor. Obviously, at NLO the probability to find a jet of high rapidity is significantly larger 
than at LO due to the extra parton that can be emitted beyond LO. This feature is illustrated 
by the .fT-factor being well below one for low rapidities, but rising continually for larger values of 
|y jet |. The impact of scale uncertainties is quite uniform over the entire rapidity range considered, 
increasing slightly for large |y jet |. 

Figure [14] shows the invariant mass distribution of the tagging jet and the Higgs boson candi- 
date, Mj^tag- The NLO-QCD effects on this observable are somewhat smaller than on the |ymax| 
distribution, but still noticeable. In analogy to the transverse-momentum distributions discussed 
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Figure 13: Rapidity distribution of the jet with the highest value of |y jet | in ep — > ejH 
at the LHeC at LO (dashed black) and NLO (solid red) [panel (a)] and associated ET-factor 
[panel (b)] as defined in Eq. (|3.5p for different values of the factorization and renormalization 
scales, /ur = /Uf = £ Q- 



0.075 - 



> 

CD 

C5 



0.050 



ST 
x 

T3 

\ 

b 



0.025 




1.05 



t 1 1 1 1 1 1 1 1 1 1 1 r 



(b) 



1.00 - 



Sh 

o 

o 

& 0.95 

I 



0.90 



f = 2 



^ = 0.5 



0.000 



200 



400 



600 



M 



H.tag 



[GeV] 



200 



400 



600 



M 



H.tag 



[GeV] 



Figure 14: Invariant mass distribution of the tagging jet and the Higgs boson candidate in 
ep — > ejH at the LHeC at LO (dashed black) and NLO (solid red) [panel (a)] and associ- 
ated .fT-factor [panel (b)] as defined in Eq. (|3,5p for different values of the factorization and 
renormalization scales, /iR = fip = £ Q. 
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above, this observable is found to exhibit scale uncertainties increasing with the kinematic in- 
variant, but still rather mild over the entire range in Mn.tag considered. 

4 Conclusions 

In this work, we have presented an NLO-QCD calculation of CC and NC Higgs production at an 
electron-proton collider. We have developed a flexible Monte Carlo program keeping full track of 
the kinematic properties of all particles involved in the reactions. We have studied the theoretical 
uncertainties being associated with the NLO-QCD predictions for inclusive cross sections as well 
as experimentally accessible observables within selection cuts for an optimized extraction of 
the Higgs signal from various backgrounds. The NLO-QCD corrections to the cross sections 
are small, changing the LO results by less than one percent for inclusive settings and about 
3% after optimized selection cuts have been applied. The relative size of the QCD corrections 
decreases when the lepton-beam energy is increased. At NLO, the residual scale uncertainties 
are at the permille level, indicating that the perturbative calculation is under excellent control. 
These findings support the conclusion of Han and Mellado [15] , based on LO studies, that Higgs 
production in ep collisions should allow for an extraction of the bottom Yukawa couplings at the 
LHeC with backgrounds and theoretical uncertainties being well under control. 
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